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ABSTRACT 
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There are a number of sources of scattering in binary optics: etch depth errors, line edge errors, 
quantization errors, roughness, and the binary approximation to the ideal surface. These 
sources of scattering can be systematic (deterministic) or random. In this paper scattering 
formulas for both systematic and random errors are derived using Fourier optics. These 
formulas can be used to explain the results of scattering measurements and computer 
simulations. 


1.0 INTRODUCTION 

For a conventional optical element, such as a mirror, scattering comes from surface roughness or 
inhomogeneities of surface constants. These surface variations are usually random and are best treated by 
statistical analysis. In some cases the manufacturing method leads to periodic or quasi-periodic surface features. 
Single-point diamond turning is an example of a manufacturing process that leaves the surface with accidental 
periodic surface heights. In these cases the surface bears some resemblance to a diffraction grating and there 
will be distinct diffraction orders. A randomly rough surface can be regarded as the superposition of an infinite 
number of diffraction gratings, each with a different grating spacing and hence a different first order diffraction 
angle. If the rms roughness of the surface is much smaller than the wavelength of light, then only the first 
diffraction order is of any consequence. 

In a binary optic the second and higher diffraction orders are important because the depth of the pattern 
is about one wavelength. In some cases these higher orders of diffraction are desirable, but in many cases 
anything except the first order can be considered "scattered" light. Scattering from binary optics will also come 
from the inevitable random surface roughness, just as it does from conventional optics. Binary optics calls for 
extremely precise positioning and control of etching to insure that each groove is the proper width and depth, 
errors of fabrication is another source of scattering. Errors of fabrication might be either systematic 
(deterministic) or random. We shall first consider systematic errors in section 2, and then random errors in 
section 3. 


2.0 SYSTEMATIC ERRORS 

Systematic errors, or deterministic errors, include: etch depth errors, line edge errors, mask alignment 
errors, quantization errors, and the shortcomings of an A/-lcvel binary optic approximation to the ideal 
continuous phase profile. We shall compare computer simulations to analytical formulas. The computer 
simulations were performed on a personal computer using the software tool Mathcad 2.5. Where possible these 
results are compared to the work of others. 

2.1 GRATINGS 

Several types of systematic errors have been investigated in the past in connection with diffraction 
gratings. When a groove position in the grating varies from the correct position, some energy will be scattered. 
Periodic variations result in the concentration of scattering into spurious diffraction orders called ghosts. Ghosts 
are common in ruled gratings because of periodic errors in turning screws and so forth. An approximate value 
of the intensity of such ghosts relative to the main diffracted order was given by Rowland in 1893. As presented 
by Hutley, 1 
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( 1 ) 


/ ghost/ / (" , ) = (*™M 2 

where m is the diffraction order, e is the maximum departure from the correct position, and a is the period of 
the grating. The location of the ghost will depend on the period of the error. As an example, if e = 0.1/Lim, a 
- l.OjLim, and m = 1, then the ghost intensity is approximately ten percent of the main diffraction order. We 
note in equation (1) that the ghost intensity is proportional to the square of the diffraction order. There should 
be no ghosts around the zero or central order. Random variations in the position of the grooves cause "grass," 
a broad spectrum of ghosts of random amplitudes in the plane of the diffraction orders. Scattering from surface 
roughness, on the other hand, has no preferred orientation and is just as intense out of the plane of incidence 
as in the plane of incidence. 

The diffraction efficiency of a grating is the fraction of the incident energy that is diffracted into a 
particular order. If only one order is desired, then the diffraction efficiency is a measure of the total unscattered 
light. We consider a diffraction grating that consists of M steps per period. All step widths and heights are 
assumed equal. We use Fourier optics theory and assume that the grating pattern extends to infinity. With these 
assumptions it can be shown that the relative diffraction efficiency 7)(p, M) of order p and with M levels is 

7)(p y M) = {[sin(7T/?/A/) sin(7Tp - <p/2)]/[ltp sin(7Tp/M - 0/2A/)]} 2 (2) 

where 0/M is the change in phase in each step. 2 We note that the highest efficiency for order p occurs when 
0/2 = 7 Tp. Equation (2) is often used for binary optics in general. The assumption is that in a small region of 
a binary optic the pattern is an approximation to a grating. 

2*2 ETCH DEPTH ERRORS 

The relative loss in diffraction efficiency with etch depth error can be obtained from equation (2). Let 
r? 0 t (P’M) th e diffraction efficiency optimized for an order p and with M steps. Let 6 be the relative etch 
depth error, then the phase is given by 

0 = 2np(l + 6). 

The diffraction efficiency is therefore given by 

T)(p,M) = {[sin(/>e7r)]/[A/ sin(p£T/A/)]} 2 . (3) 

Figure (1) is a plot of the diffraction efficiency for Fresnel zone plate design FZP3A as a function of 
percent error in etch depth. Table 1 summarizes the design parameters for several Fresnel zone plates. The 


TABLE 1. FRESNEL ZONE PLATE DESIGNS 


Design 


Levels 

T 

/(mm) 

T# 

n 

b( mm) 

FZPl 

0.6328 

16 

316 

10 

10 

1.457 

0.35 

FZP2A 

0.6328 

4 

79 

10 

10 

1.457 

0.35 

FZP2B 

0.6328 

4 

2005 

254 

10 

1.457 

00 

FZP3A 

0.6328 

2 

40 

10 

10 

1.457 

00 

FZP3B 

0.6328 

2 

40 

10 

10 

1.457 

0.35 

FZP3C 

0.6328 

2 

395 

100 

10 

1.457 

3.50 

FZP4 

10.6 

16 

19 

10 

10 

4.0 

0.35 

FZP5 

10 

16 

3360 

152.4 

3 

4.0 

00 

FZP6 

10 

16 

300 

12.7 

1 

4.0 

00 

FZP7 

10 

8 

1680 

152.4 

3 

4.0 

00 

FZP8 

10 

8 

300 

12.7 

1 

4.0 

00 


data points were calculated using the approximate Rayleigh-Sommerfeld equation which for a circularly 
symmetric Fresnel zone plate, on axis, can be written as 
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where, 


U(0,0, z) = s. (2nz/\i) 0j / aJ ai + 1 r/ir 2 + z 2 ) exp[/7c(z 2 + r 2 )*] dr 


( 4 ) 


= [2jj\/M + (A;/Af)¥ ( 5 ) 

0^ = exp{-27T/[l + e]mod(/, M)/M} (6) 

and z is the distance from zone plate to observation point, / is the focal length, T is the maximum number of 
transition points A/ is the number of levels, and e is the fractional etch depth error. We define the function 
mod(/', M) as the remainder when j is divided by M. For this example the illuminating beam was assumed to be 
uniform, so the Gaussian beam radius b (at the 1/ e 2 points of intensity) is infinite. 

For an etch depth error that is uniform across the surface we can use equation (3) for the far-field, 
paraxial case. In a 2-level Fresnel zone plate an etch depth error of e will alter the phase over half the surface 
of the optic (in the grooves only), so we use e/2 for the error term in equation (3). This curve, also shown in 
figure (1), matches the data very well. Data similar to figure (1) was shown by Cox, et. al. 

Uniform etch depth error is one type of systematic error that can occur in a binary optic. The etch 
depth can also vary across the binary optic; for example each etch depth may be 5 percent deeper than it should 
be near the center, gradually becoming 5 percent shallower near the edges. This might occur when the narrower 
grooves near the edges reduce the etching rate. We can model this type of error by writing 

<p . = exp{-27Ti[l + e(l - 2tfj/fl T )]mod(/, M)/M) (7) 

for the phase term 0. in equation (4). In computer simulations using the Huygens-Fresnel integral (the 
approximate Rayleigh-Sommerfeld equation), Goodman and Farn found empirically that the diffraction 
efficiency of a Fresnel zone plate falls off as the square of the etch error. The effect of etch errors is only a 
weak function of the number of masks, and it is independent of the f-number and radius. They considered 
Fresnel zone plates with f-numbers // 1 and // 3 and radii of one and three inches. 

Figure (2) shows results identical to Farn and Goodman's data for a 16-level and an 8-level Fresnel 
zone plate. The relative diffraction efficiency is the peak intensity (equations (4), (5) and (7)) divided by the 
product of the peak intensity of a perfect lens and the maximum efficiency of an M-level zone plate (equation 
(2)). Designs FZP5, FZP6, FZP7, and FZP8 from Table 1 are used in the calculations. The results do not 
depend on the f-number, and depend only slightly on the number of levels (hence on the number of masks in 
the fabrication process). Also shown in figure (2) are curves calculated from equation (3), where £ = 7/8 e for 
an 8-level zone plate, and £ = 15/16 e for a 16-level zone plate. The relative etch depth error £ appears to be 
the average systematic etch depth error along an axis of the Fresnel zone plate. 

A systematic etch depth error decreasing from center to edge should shift the focal position. When 
the efficiency is calculated at the optimum focal position, then the diffraction efficiency decrease is not nearly 
as severe as calculated by Farn and Goodman in reference (4). The upper data points in figure (2) shows my 
calculations for diffraction efficiency at best focus (empirically determined). A good fit to the data is obtained 
by using equation (3) with £ = ejj\ which is the mis value of the systematic etch depth error e along an axis 
of the Fresnel zone plate. At best focus it is found that the diffraction efficiency is essentially independent of 
the number of masks, the f-number, and the radius. 

23 LINE EDGE ERRORS 

One systematic line edge error would be a linear increase in the position of the transition points. We 
write the focal length / from equation (5) in terms of the ;th transition point flj, the wavelength X, and the 
number of etch levels M, 

f = a 2 M/2jX - j\/2M. (8) 
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A uniform expansion of the transition points {#j} by the factor B results in a new focal length / ' given by 

/' = B 2 a 2 M/2jX - ;\/2M 

= B 2 f + (B 2 - l)Xj/2M. (9) 

In the paraxial case (i.e., when j is small) the second term of equation (9) can be ignored and the new focal 
length is just B 2 times the old focal length. In general the focal length depends on j , which means there is 
spherical aberration. The scattering in this case Is a spreading of the energy around each diffraction order; and 
the larger the aperture, the larger the spreading: 

Over- or under-exposure of the photoresist during the pattern transfer process can result in systematic 
line edge errors. Under-exposure results in the grooves being narrower than they should be. Farn and 
Goodman found that diffraction efficiency (defined as the Strehl ratio) of a Fresnel zone plate falls off linearly 
with the line edge error, the effect of line edge error is a strong function of the f-number, a weak function of the 
number of masks and Is independent of radius. 4 

Cox et. al. measured the diffraction efficiency of an// 10 Fresnel zone plate in ten "zones" of equal 
area. 3 The first zone is the central portion of the zone plate, covering ten percent of the entire zone plate. Two 
zone plates were deliberately under-exposed so that the etched portions were 1.35 /im narrower than they should 
be. A plot of the measured local diffraction efficiency (fraction of power) is shown in figure (3) for the average 
of the two zone plates. Also shown is a plot of the calculated local diffraction efficiency in each of the ten 
"zones" as defined by Cox, et. al. 

We expect the diffraction efficiency to decrease with increasing zone number since the error e becomes 
relatively more significant as the spacing decreases. In a 2-level zone plate a line edge error represents an area 
of the zone plate that does not contribute at all to the peak intensity. We therefore subtract the error fraction 
of a period from the efficiency. In a 2-level zone plate a period is every two transition points. In design FZP3A 
there are forty transition points, hence four transition points in each of ten "zones" as defined by Cox, et. al. If 
(I' is the average period in the ith "zone", or 

d \ = K+2 * °4i + °4i + 4 ’ °4i + 2)/ 2 ’ ( 10 ) 

then the local diffraction efficiency t] ■ is approximately given by 

r).{e) = (1 - e/d.) [M/ir sm(n /M)f, (11) 

where e is the line edge error. When e = 0 equation (11) reduces to the standard formula for diffraction 
efficiency with no errors (equation (2)). To compare with the Cox data we use the value of e = 1.35 /xm 
measured by Cox, and we find a good fit using equation (11), as shown in figure (3). The diffraction efficiency 
is seen in equation (11) to decrease linearly with line edge error and is a strong function of the f-number (the 
smaller the f-number, the smaller the value of d in the outer zones) as observed by Farn and Goodman. 4 

2.4 MASK ALIGNMENT ERRORS 

The experimental and theoretical work of Cox 3,5 shows that a fixed lateral error, or a mask alignment 
error, leads to a fall off in diffraction efficiency that is more severe in the outer zones of the zone plate than near 
the center. In the outer zones a fixed error represents a greater percentage of the pattern and we expect it to 
cause the greater efficiency decrease. By measuring the local diffraction efficiency of a number of zone plates 
with design FZP2B, Cox found that the fall off in efficiency is approximately linear from zone to zone. I have 
plotted in figure (4) the diffraction efficiency in the outer most (i.e. 10th) zone as a function of the mask 
alignment error. The data points with a small alignment error represent typical fabrication errors, the large 
alignment error values were intentional. The straight line in figure (4) is given by 

n, 0 (e) = (1 - e/Aa 10 ,ave) t A/ / 7r sin(7T/A/)] 2 (12) 
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where M = 4 for this 4-level design, and Aa 10 = 3.25 Hm is the mean separation of transition points in the 
tenth zone. The mean separation is much smaller than the average period used in equation (11). Cox concludes 
that mask alignment error is more significant than etch depth error or line edge error. 

2.5 M-LEVEL APPROXIMATION ERRORS 

A 100 percent efficient Fresnel zone plate would have a continuously curving profile. Equation (2) 
indicates that a 16-level binary optic approaches that level of efficiency. How is the energy distributed that is 
not diffracted into the first order? We know that in the absence of fabrication errors the energy goes into other 
diffraction orders which will have very low intensity at the first-order focal plane. To calculate the scattering we 
first reduce the diffraction "ringing," caused by the finite size of the Fresnel zone plate. The amount of energy 
outside of the central lobe can be minimized by assuming that the Fresnel zone plate is illuminated by a 
Gaussian beam that has low intensity at the edges of the aperture. 

Figure (5) compares the intensity as a function of angle for a cylindrical Fresnel zone plates with 2-, 
4- and 16-levcls. We are using designs FZP3B, FZP2A, and FZP1. The intensity is relative to the peak intensity 
/ of a perfect lens with the same f-number. The data has been smoothed by averaging over sixteen successive 
data points. The first point plotted is the mean of the first sixteen computer generated data points, I. to / 16 , 
divided by I . The second point plotted is the mean of / 2 to / 17 divided by / p , and so forth. There are 2048 data 
points for each curve, and the points are spaced at intervals of A. 

We notice in figure (5) that the intensity of the scattered light from the 2-level zone plate seems to 
drop to the intensity of the 4-level zone plate at about 5.7 degrees. This can be explained as follows. In the 1st- 
ordcr focal plane the biggest visible difference between the 2-level scattering and the 4-level scattering will be 
the contribution of the 3rd-order energy. The diffraction efficiency equation (2) indicates that for a 2-level 
binary optic about 4.5 percent goes into the 3rd-order. The 4-level (and higher) binary optic has no energy in 
the 3rd-order. 

The paraxial focus of the 3rd-order is at//3, where / is the lst-order focal length. Using geometrical 
optics the extreme ray from the edge of the zone plate which passes through the 3rd-order focus strikes the 1st- 
order focal plane at x Q ' given by 

x 0 7(f-// 3 ) = (j/2t*)/m) ( 13 > 

where 1* is the lst-order f-number of the zone plate. Hence 

v = //I*. < 14 ) 

This is the geometric edge of the 3rd-order diffraction in the lst-order focal plane. As seen from the center of 
the zone plate this edge makes an angle O', given by 

tan = x Q '/f 

= 1/1*. (15) 

For an//10 zone plate O' = 5.7 degrees, in good agreement with our observation in figure (5). 

3.0 RANDOM ERRORS 

Random errors are processing errors that are statistical in nature. In this section we discuss some 
random errors that have been noted for gratings. We then derive a formula for random roughness for binary 
optics. We show that this formula can also be used to calculate the scattering from random etch depth errors 
and random line edge errors. We compare the formula to computer simulations of random fabrication errors. 
We use Dammann gratings and Fresnel zone plates as examples. 
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3.1 GRATINGS 

From the theory of gratings a random change in form from groove to groove is called "accidental error 
of amplitude," and is a form of random error independent of grating order. 6 A random error in the position 
of the groove causes "grass," as mentioned previously. The proportion of the energy scattered into grass is given 
by 1 

W/ 1 ^) = ( 47r " ,e rms/^ 2 ( 16 ) 

where e fms is the mis error in groove position. According to this formula grass increases as the square of the 
diffraction order, the same as the intensity for ghosts. 

32 RANDOM ROUGHNESS 

We consider now the problem of a binary optic with slightly rough surfaces. 7 We shall consider the 
binary optic to behave like a grating and use the Fraunhofer diffraction formula. The illumination of the grating 
will be by a collimated laser beam of wavelength A and beam radius b. The beam waist is at the grating. The 
field at the point (x Q , y 0 , z ) is given by (ignoring some unimportant phase terms) 

U(v 0 > v 0 ’ z) = 1/A z f JD a f J0 '° exp (lirimxja) - 1)] 

exp[-/A:(r 0 r 1 + ytf^/z) expK*, 2 + y 2 )/b 2 ] dx y dy , (17) 

where E/4 m exp(2? Timx^/a) is the Fourier series of the field induced by the binary optic in the plane, 

a is the grating period, n is the index of refraction of the grating material, and h is the roughness or surface 
height variation function. We multiply U(r Q , y Q , z) by its complex conjugate to obtain the intensity. 

We next assume that the distribution of surface heights h is random with a Gaussian probability; 
therefore 8 

<exp{l/[h(jr 1 ,y 1 ) - h(* 2 , y 2 )](/i - 1)}> = exp{-[fc(/i - 1 )ct] 2 [1 - R(r, - x 2 ,y, -y 2 )]} (18) 

where < > denotes the average or expected value, G is the mis surface roughness and RQ is the autocorrelation 
function. If we assume that the surface roughness is small, ko < < 1, then 

exp{-[^(/i - l)a] 2 [1 - R(r 1 - x 2 ,y, -y 2 )]} ~ 1 - [k(n - l)a] 2 [1 - R(r, -x 2 ,y y -y 2 )]. (19) 

Next, assume that the autocorrelation function is exponential, 

R (*i - x 2< y\ -y^} = ex p{-[(^i ' x 2> 2 + (y-i -yz> 2 \ y2 / j } ( 20 ) 

where T is the correlation length. Recognizing that we now have the zero order Hankel Transform of the 
autocorrelation function we find that 

< I(-*c - v 0’ z ) > = < ‘rough*** - V 0’ z )> + < ‘smooth ( x 0’ ^0’ z ) > ( 21 ) 

where 

< 'rough ^Wo’ 2 ) > = l naTbk ( n ' !)/^] 2 S M m | 2 I 1 + [(2 nnt/a - /zf + (y^/z) 2 \T 2 y 3/2 , (22) 

<I smooth( x 0’>'0’ 2 ) > = (b 2 n/yzf {1 - [*(« - l)cr) 2 } £M m | 2 exp{-2Z> 2 jr 2 f(v 0 /Az - »i/a) 2 + (y 0 /Az) 2 ]}. (23) 

Wc can compare these results to the scattering formulas of others by rewriting this expression as the 
bidirectional transmission distribution function (BTDF). The sine of the specular angle & 2 m * n *^e order 
is given by 
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(24) 


sin0 2m = Xrn/a 

and the sine of the scatter angle is, to the approximations we are using, 

sin0 s = x Q /z. (25) 

We define the spatial frequency p m by 

Pm = *(sin* s - sin(? 2 m ). (26) 

Along the axis, y Q = 0, and we can write the BTDF using equation (22), and equations (24) to (26) by 

BTDF = (2* 4 /7T) [OT(n - l)/2] 2 cos 0, Z \AJ 2 /[I + P m 2 T 2 f /2 . (27) 

The typical two-dimensional bidirectional reflection distribution function (BRDF) with an exponential 
autocorrelation function, for normal incidence is given by® 

BRDF = (2k 4 /n) (OT) 2 cos 9 t F(0 $ ) /[I + p m 2 T 2 ] 3/2 (28) 

where F((? ) is approximately one for small angles, and is exactly one in some theories. Comparing equation (28) 
with equatton (27) we see that one difference is the factor [(/i - l)/2] 2 . This difference exists because the phase 
change on reflection from surface features is proportional to twice the height of the feature, but in transmission 
the phase change is (n - 1). The second difference is the term H\A m \ 2 . Evidently the existence of the grating 
replicates the scattering pattern around each diffraction order. The amplitude of the scattering is scaled by the 
amplitude of the diffraction order. 

In a similar manner we can derive a formula for one-dimensional random roughness. We shall show 
that this formula can also predict the scattering from Dammann gratings with random etch depth and random 
line edge errors. In addition the formula will predict scattering for cylindrical Fresnel zone plates with random 
etch depth or line edge errors; the results agree well with the average of computer simulations. Equation (22) 
will predict scattering from circular Fresnel zone plates with random etch depth or line edge errors. 

The one-dimensional random roughness formula is derived by assuming that the equations are 
independent of the y 1 -axis. Instead of a zero order Hankel Transform of the two-dimensional autocorrelation 
function, we have the Fourier Transform of the one-dimensional autocorrelation function. The end result is 

< K^o’ l'o’ > = < ^rough^O’ 1'0’ ^ ^ + < ^smooth^O’ ^0’ > (^) 

where 

^roughCWo- 2 ) > = W" • l)o/Xz\ 2 brJ(nr)T.\A m \ 2 [\ + T 2 (kx 0 /z - 2nm/a)Y (30) 

< Smooth (WO’ z )> = (V^) 2j r{l - Wn - l)af}i:\A m \ 2 exp{-2b 2 n 2 \x 0 /Xz - m/a] 2 }. (31) 

We recall that O is the nns roughness, T is the correlation length, k = 27 T/X, b is the radius of the Gaussian 
beam, n is the refractive index of the substrate, a is the period of the grating, and M m | 2 ( as lt turns out ) is the 
diffraction efficiency of the /nth order. 

33 ETCH DEPTH ERRORS 

We now consider etch depth errors where the etch depth varies randomly from groove to groove. We 
can consider Dammann gratings and cylindrical Fresnel zone plates as one-dimensional rough surfaces. Circular 
Fresnel zone plates are two-dimensional. 
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33.1 DAMMANN GRATING 

Dammann gratings are convenient to consider because they are simple two-level binary optics where the 
diffraction orders are spatially separated. We consider a number of designs, as shown in Table 2. 


TABLE 2. DAMMANN GRATING DESIGNS 


Design 

Transition points 

depth 

period 

A 

Equal intensity 


° 1 (M m ) o 2 (Hm) 

/i(jim) 

a (Jim) 

(Mm) 

(orders) 

DGl 

43.02 86.96 

0.545 

200 

0.633 

0, ±1, ±2, ±3 

DG2 

43.02 86.96 

9.129 

200 

10.6 

0, ±1, ±2, ±3 

DG3 

10.75 21.74 

0.545 

50 

0.633 

0, ±1, ±2, ±3 

DG4 

13.74 30.06 

1.23 

60.12 

10.6 



Figure (6) shows the relative intensity of Dammann grating design DG1 with a random etch depth 
error. In addition to the values in Table 2, b = 4.5 mm, / = 1.82 m, and n = 1.45702, except for DG4 where 
b = 2.5 mm and n = 3.4178. The standard deviation of the error is S =~UJB 1 Jim . The intensity profile was 
calculated with four sets of random numbers and the results were averaged together. The results were also 
smoothed by averaging over sixteen data points, and normalized by dividing by the peak intensity in the error- 
free case. In Figure (6) the zero diffraction order is afzerd degrees^ and the next four posltive diffraction orders 
are also shown. The negative diffraction orders would be similar and are not plotted. 

Also shown in Figure (6) is a plot oT equation (29). file correlation length T we estimate to be 
approximately equal to the width of the etched portion of the Dammann grating. In this portion the correlation 
is unity, falling off to zero outside of the etched area. For design DG1 there are two etched portions per period; 
the width of each is 86.96 - 43.02 = 43.94, so we set T = 44. The mis roughness O is 

a = (< S 2 r/a + 0 + 6 2 T/a + 0)* 

= s (IT /of. (32) 

In Figure (7) the standard deviation of the etch error has been increased to 6 = 0,1 Jim. The average 
of four computer simulations using the Fresnel diffraction formula is shown together with the 1-D roughness 
equation (29). The value of T is again 44 Jim, and we use equation (31) to calculate the mis roughness. It 
appears to be a good fit. We see from equation (29), and it is conFirmed in the computer simulations illustrated 
in Figures (6) and (7), that the scattering is proportional to <5 2 , since S and O are linearly related through 
equation (31). The diffraction efficiency can be given by equation (31), which indicates that for random etch 
depth errors with two etched surfaces per period, the diffraction efficiency decreases by 

n - V 1 * W" - 

= r] Q {\-(2T/a)[k(n-l)&\ 2 } (33) 

using equation (32). The efficiency with no errors is rj Q , and is given by equation (2). 

What effect does wavelength have on the scattering? Consider design DG2, which has the same 

transition points (hence the same pattern), but differs in the design wavelength and etch depths The relat ive 
intensity profile is shown in figure (8), along with a plot of the one-dimensional scattering formula (29). Formula 
(29) agrees as well with computer simulations for the 10.6 Jim wavelength as for the 0.6328 Jim wavelength. The 
etch depth error has a standard deviation of 0.1 Jim, the same as the error illustrated in figure (7), but the 
relative scattering is much less. This is in keeping with the prediction of equations (30) and (31) that the relative 
scattering is 

< W ( Wo’ 2 ) > / I peak = ' 2)0 /\zfbT J{lTT) X\AJ Z [ 1 + T 2 (kx Q /z - lllm/afx' /{b/\zfn\A^\ 2 
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= [k(n - l)o] 2 Tj2I,\A m /A peak \ 2 [bJn + bjm 2 (kxjz - 27 rm/a) 2 ]' 1 (34) 

assuming there is not much overlap of energy between diffraction orders. We see from equation (34) that the 
relative scattering is proportional to ( ko ) 2 and therefore proportional to (6 /A) , as we can also see by comparing 
figures (8) and (7). 

3 32 CYLINDRICAL FRESNEL ZONE PLATES 

Continuing with our examination of random etch depth errors, we now look at cylindrical Fresnel zone 
plates. This type of zone plate will show one-dimensional scattering characteristics. The cylindrical Fresnel 
zone plate creates a line focus (at Xq — 0) for each diffraction order at a distance from the zone plate that 
depends on the diffraction order. We shall concern ourselves only with the scattering from roughness in the first 
diffraction order; in the lst-order focal plane scattering from other diffraction orders is insignificant. The relative 
scatter intensity from roughness is therefore (using equation (34)), 

< WWo* z ) > / , peak « W n ' Wfrjl [bjn + bJ*{Tkx Q /z) 2 \ A (35) 

for the cylindrical Fresnel zone plate. 

As before, the correlation length T is the mean width of the etched portions; 

T = mean{aj + 1 - flj} (36) 

which for design FZP1 gives T = 1.58 jJ.m. For the Fresnel zone plates all levels were assumed to have some 
etch depth error, so the mean is taken over all intervals. This also implies that 0=6. Figure (9) should be 
compared to figure (10) where the standard deviation of the etch depth error is 6 = 0.1 flm and 6 = 0.01 M m, 
respectively. The scattering drops two orders of magnitude with one order of magnitude decrease in the error. 
Again equation (35) provides a good prediction of the average value of the scattering at a particular angle. 
Equation (35) predicts that the scattering is proportional to CT , and o 2 = 6 2 in this case. 

A 4-levcl Fresnel zone plate (FZP2A) shows behavior similar to the 16-level zone plates just considered. 
Figure (11) illustrates this case when the standard deviation of the etch depth error is 6 = 0.1 M m . A plot of 
equation (35) is shown for comparison. Equation (36) gives T = 6.33 ^tm; and 0=6. The random scattering 
in this case is less than the systematic error for angles greater than about three degrees; averaging is necessary 
to reveal the scattering. The gaps in the curve indicate where the average intensity with etch depth errors is less 
than the intensity without etch depth errors. Further averaging should make the relative scatter intensity 
converge to the plot of equation (35). 

3 33 CIRCULAR FRESNEL ZONE PLATES 

We now consider random etch depth errors in circular Fresnel zone plates. In this case we must use 
the two-dimensional roughness scattering formula, equation (21). For our purposes here we want the scatter 
intensity from the "rough" part divided by the peak intensity of the smooth part in the absence of roughness. 
Using equations (22) and (23) we find 

< Wwo-^/W = w* - ] ) aT M 2 1 1 + ( Tfc * 0 / 2 ) 2 + ( t *V 2 ) 2 ]‘ 3/2 ( 37 > 

for the two-dimensional or circular Fresnel zone plate. The mis roughness O will be, as before, equal to the 
standard deviation of the etch depth error, S. The correlation length will again be given by equation (36). 

The relative scattering for the circular Fresnel zone plate with design FZPl is shown in figure (12). 
The relative scattering in the two-dimensional case is much less than the relative scattering in the one- 
dimensional case (figure (9)) for the same average etch depth error. 

Design FZPl also has much less relative scattering than design FZP4, shown in figure (13). What is the 
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reason for this? There are three terms in equation (37) which change from design FZP1 to FZP4; T, A, and 
n. The correlation length T in turn is a Function of the wavelength A, the number of levels M y and the f-number 
f*. We now derive an approximate value for T. If the f-number is not too small we can approximate the 
transition point a j from equation (5) by using the paraxial approximation, 

Oj 2 * 

Hence, 

°j+i 2 ■ a \ ~ zx f/ M 

and therefore 


a j + l - « 2A//(Afa j + 1 +Map. 

We have then that 

r = mean(aj + 1 - a p ~ [2A//A/] mean(l/[jj + 1 + a.]) = 2\f/Ma T 

r « 4\f/M. (38) 

We see that the term /c r in equation (37) is approximately the same in both design FZP1 and FZP4, hence the 
difference in relative scattering levels is due solely to the difference in retractive index n of the two designs. 

3.4 LINE EDGE ERRORS 

We now examine the effects of random line edge error in Dammann gratings and Fresnel zone plates. 
The location of the line edge is made random by adding e- to the transition point where e- is a random 
number with mean of zero and standard deviation <5. Our random roughness formulas were derived under the 
assumption of a small random phase shift to the surface, but we will attempt to use these formulas for the 
random line edge errors. 

A random line edge error could be regarded as a perturbation from the ideal pattern. The height of 
the perturbation is the step height, h. The mis roughness O is then the square root of the average amount of 
height perturbation, or 

O = \m<w>h 2 /a J 1 / 2 (39) 

where nx is the number of edges which vary from their proper positions in a period a , and <w> is the mean 

width of the perturbations. For a uniform probability over either the interval [0, 6j 3] or [-&J 3, 0], 

<h» = &(/. 3)/2 (40) 

and & is the standard deviation of the error over that interval. 

In the derivation of the mis roughness formula we assumed a two-sided exponential probability density 
function, (1/2t) exp(- \x\/t). Hence the standard deviation S is related to the correlation length r by 

<5 = [(1/2T) J^Vexp(-M /T)dx)^ 2 

.. . = tjz _ m 

With a and T now given by equations (39), and (41), we can compare the computer simulations of random line 
edge error to plots using the roughness formulas. 
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3.4.1 DAMMANN GRATINGS 

We begin our examples of random line edge errors with the Dammann gratings. Figure (14) shows the 
relative intensity versus angle for Dammann grating DG1 with line edge error S - 1 fi m. A plot of the 
roughness formula (29) is also shown (it is almost constant at 1(T). There are four edges per period in design 
DG1, so m = 4 in equation (39). The step height h is given in Table 2. The scattering predicted by equation 
(29) appears to be a little high compared to the computer simulation, possibly because the assumption made in 
the derivation that h is small compared to A is not very good. Figure (15) shows the case for 6 = 0.1 Jim. 
Equation (30) correctly predicts an increase in scattering of magnitude 100 when 6 is increased by a factor ot 
10. To see this in equation (30) we note that 

<I rough( jr 0^0’ 2 ) > 

<x <w> r 

« S 2 (42) 

using equations (39), (40), and (41). When the period is decreased by a factor of four (design DG3), the 
scattering increases by a factor of four (not shown here). The scattering is proportional to O in equation (30), 
and therefore proportional to 1/a, according to equation (39). 

When the wavelength is increased, and the step height is increased proportionally (design DG2), the 
relative scattering intensity remains the same for random line edge errors. The roughness formula (34) predicts 
this since the relative scattering is proportional to (ok) 2 , which is proportional to (/i/A) by equation (39). 

Figure 16 is a comparison of the roughness formula with measured data. The design is the same as 
DG1 except the etch depth was (accidentally) 0.662 (ini instead of 0.545 flm. This results in a reduction of the 
zero order about a factor of ten from the design value. Measurements using a WYKO profilometer indicate a 
mis roughness of about 0.43 nm, too little to explain the measured scattering. The random etch depth error was 
measured to be about 1.0 nm, which could account for some of the scattering. A random line edge error with 
standard deviation of 0.025 (tm, when added to the etch depth error, gives a good approximation to the measured 
scattering. This level of random line edge error was too small to be measured with the WYKO, and is less than 
a typical manufacturing error. 

Another comparison of measured data and the roughness formula is shown in figure (17) for design 
DG4. The mis roughness of the lands was about 1.0 nm, but the grooves had a slope giving an effective 
roughness of 27.6 nm. Since this slope was periodic we would expect a systematic error, rather than a random 
error. Plotting these values of etch depth roughness shows 1.0 nm accounts for only a small part of the 
measured scattering, and 27.6 nm predicts orders of magnitude more scattering than was measured. A line edge 
error with standard deviation of 0.03 jitm agrees reasonably well with the scattering measurement. 

3.4.2 CYLINDRICAL FRESNEL ZONE PLATES 

We now consider random line edge errors in Fresnel zone plates. Figure (18) is a plot of the relative 
scatter intensity of a cylindrical zone plate, design FZP1. The systematic error part of a 16-level zone plate has 
been subtracted out. We add a random number e ■ to each transition point a-. To plot the roughness equation 
we use the same formulas for T and a as we did for Dammann gratings. We interpret equation (39) using m 
= T for the number of transition points or edges in the half-width or radius a of the design. We could write 

O = h\6(j3)T{#/f\ h (43) 

using equations (39) and (40). The relative scattering increases by a factor of 100 when the line edge errors 
increase by a factor of 10, as it does for Dammann gratings. 

Design FZP2 has one fourth the transition points of design FZP1, but the step heights are four times 
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as high, since 


h = \/[(n - 1 )M\. 


(44) 


Roughness is proportional to hjT \ as we see from equation (43), so design FZP2 has twice the roughness of 
design FZP1, and therefore four times the scattering for the same random line edge error. Computer 
simulations confirm this prediction of the roughness formula. 


There is less scattering from design FZP4 than design FZP1. It appears to be a result of differences 
in the factors A, T y and n. The number of transition points T is reduced approximately in proportion to A. The 
step height is increased in proportion to A, but decreased in proportion to (n - 1). We have from equations (30), 
(43), and (44) 

tough « T(« - 1)CT/A] 2 

a (/j - l) 2 h 2 T/X 2 
« (n - l) 2 [\/{n - l)] 2 1/A 3 

« 1/A. (45) 

We see that the refractive index n doesn't really affect the scattering intensity for line edge errors, and the net 
result is that the scattering is inversely proportional to the wavelength. However, if one particular design were 
used at two wavelengths, then T and o would be constant, and the scattering would be proportional to 1/A 2 . 

3.4.3 CIRCULAR FRESNEL ZONE PLATES 

Lastly we consider designs FZPl and FZP4 for random line edge errors in circular Fresnel zone plates. 
The intensity is found using random numbers e- added to the transition points a-. These results are compared 
to the roughness formula (37). As before, the correlation length is <5//2, and the mis roughness G is the mis 
surface height with respect to the design step height. For a circular Fresnel zone plate 

a = {SjKa, + <w>) 2 - flj 2 ](/ 1 /a T ) 2 } 1 /2 
= (h/a T )[T<w> 2 + 2<w>E j a j ] 1/2 

= (h/a r )[T6 2 3/4 + <573 S^] 1 / 2 . (46) 

Figure (19) shows that the relative intensity of the scatter is down six or seven orders of magnitude 
from the peak intensity, even for a fairly large line edge error of 1 [Jim standard deviation. As in the one- 
dimensional case, the roughness formula predicts a decrease in relative scattering proportional to 1/A for 
optimized designs. 

4.0 CONCLUSIONS 


Using Fourier optics methods some formulas for scattering from binary optics have been developed. 
These formulas are similar to standard random roughness formulas used to predict scattering from smooth 
surfaces. The formulas are shown to agree reasonably well with computer simulations and actual measurements. 
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Etch depth error (%) 



FIGURE 1. Diffraction efficiency versus uniform 
etch depth error (%), Fresnel zone plate (FZP3A). 
Computer simulation (X) compared to equation (3). 


FIGURE 2. Relative diffraction efficiency versus 
systematic error (%). 16-lcvel (0), 8-level (X), 
and best focus ( + ); compared to equation (3). 
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FIGURE 3. Diffraction efficiency versus zone 
number. Comparing Rayleigh-Sommerfeld (X), 
equation (11), and data (-0-0-) from ref. (3). 



FIGURE 4. Diffraction efficiency in the outer zone 
versus mask alignment error (fJim). Equation (12) 
compared to measured data (X) from ref. (5). 




FIGURE 5. Relative scatter intensity versus 
scattering angle. Cylindrical Fresnel zone plate; 
16-level (bottom), 4-level, and 2-level (top). 


FIGURE 6. Relative intensity, Dammann grating 
(DG1). Etch depth error S = 0.01 jLtm. Compared 
to 1-D roughness scattering formula. 
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FIGURE 7. Relative intensity, Dammann grating 
(DG1). Etch depth error S = 0.1 [im. Compared 
to 1-D roughness scattering formula. 



FIGURE 9. Relative scatter intensity, cylindrical 
Fresnel zone plate (FZP1). Etch depth error S = 
0.1 /xm. Compared to 1-D scattering formula. 


FIGURE 8. Relative intensity, Dammann grating 
(DG2). Etch depth error <5 = 0.1 fJ,m. Compared 
to 1-D roughness scattering formula. 



FIGURE 10. Relative scatter intensity, cylindrical 
Fresnel zone plate (FZP1). Etch depth error S = 
0.01 pm. Compared to 1-D scattering formula. 
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FIGURE 11. Relative scatter intensity, cylindrical 
Fresnel zone plate (FZP2A). Etch depth error S = 

0.1 jLlm. Compared to 1-D roughness scattering formula. 



FIGURE 12. Relative scatter intensity, circular 
Fresnel zone plate (FZP1). Etch depth error 6 
0.1 jJLm. Compared to 2-D scattering formula. 



FIGURE 13. Relative scatter intensity, circular 
Fresnel zone plate (FZP4). Etch depth error 6 = 
0.1 /Zm. Compared to 2-D scattering formula. 


132 







Angle (degrees) Angle (degrees) 

FIGURE 14. Relative intensity, Dammann grating FIGURE 15. Relative intensity, Dammann grating 

(DG1). Line edge error 6 = 1.0 Jim. Compared (DG1). Line edge error S = 0.1 /im. Compared 

to 1-D roughness scattering formula. to 1-D roughness scattering formula. 



Angle (degrees) Angle (degrees) 

FIGURE 16. Measured relative intensity; compared FIGURE 17. Measured BTDF compared to 1-D 

to 1-D scatter formula with etch depth error <5d = scatter formula with line edge error S = 30 nm. 

1 nm, and line edge error S\ = 25 nm. 
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FIGURE 18, Relative scatter intensity, cylindrical 
Fresnel zone plate (FZP1). Line edge error 6 = 1 
/Lim. Compared to 1-D roughness scatter formula. 


FIGURE 19. Relative scatter intensity, circular 
Fresnel zone plate (FZP1). Line edge error 6 = 1 
jUm. Compared to 2-D scattering formula. 
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